Dynamics of Non-Conservative Voters 
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We study a family of opinion formation models in one dimension where the propensity for a voter 
to align with its local environment depends non-linearly on the fraction of disagreeing neighbors. 
Depending on this non-linearity in the voting rule, the population may exhibit a bias toward zero 
magnetization or toward consensus, and the average magnetization is generally not conserved. We 
use a decoupling approximation to truncate the equation hierarchy for multi-point spin correlations 
and thereby derive the probability to reach a final state of j consensus as a function of the initial 
magnetization. The case when voters are influenced by more distant voters is also considered by 
investigating the Sznajd model. 
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It often happens that individuals change their atti- 
tudes, behaviors and/or morals, to conform to those of 
their acquaintances. Perhaps the simplest description of 
this conformity is the voter model [l|] , where each node of 
a graph (i.e., the social network) is occupied by a voter 
that has one of two opinions, f or |. In the voter model, 
the population evolves by: (i) picking a random voter; (ii) 
the selected voter adopts the state of a randomly-chosen 
neighbor; (hi) repeating these steps ad infinitum or un- 
til a finite system necessarily reaches consensus. Figu- 
ratively, voters have no self confidence and merely fol- 
low one of their neighbors. With this dynamics, a voter 
changes opinion with a probability pf that equals the 
fraction / of disagreeing neighbors. This proportionality 
rule leads to the conservation of the average opinion in 
the system, a feature that renders the voter model sol- 
uble in all dimensions P, However, the specific rule 
Pf = f is one among many possible and socially plausible 
relations between pj and /. 

In this work, we generalize the voter model so that 
Pf depends non-linearly on /. Non- linear voter models 
have been discussed previously, primarily by numerical 
simulations in two dimensions |3|. Here we focus on one 
dimension, where the range of possibilities for the non- 
linearity is limited. In one dimension, a voter may be 
confronted by 0, 1, or 2 disagreeing neighbors. It is nat- 
ural to impose po = 0, so that no evolution occurs when 
there is local consensus. Then the most general descrip- 
tion of the system requires two parameters, p\ and P2 
(Fig. [T|) . One parameter, which we choose to be p\ , deter- 
mines the overall time scale and is thus immaterial. The 
only relevant parameter then is 7 = Pilp\. When 7 = 2, 
one recovers the classical voter model. When 7 > 2, the 
combined effect of two neighbors is more than twice that 
of one neighbor. Equivalently, voters can be viewed as 
having a conviction for their opinion and that strong peer 
pressure is needed to change opinion. As 7 — > 00, voters 
only change opinion when are confronted by a unanimity 
of opposite-opinion voters [J]. In contrast, when 7 < 2, 
one disagreeing neighbor is more effective in triggering an 
opinion change than in the classical voter model. When 
7 = 1, one recovers the vacillating voter model 0, Q 



where voters change opinion at a fixed rate if either 1 or 
2 of their neighbors disagree with them. Finally, 7 < 1 
corresponds to a "contrarian" regime where a voter is 
less likely to change opinion as the fraction of neighbors 
in disagreement increases. 
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FIG. 1: Update illustration. A random voter changes state 
with probability p\ if it has 1 disagreeing neighbor (left), and 
with probability p2 if it has 2 disagreeing neighbors (right). 

We term 7 the conviction parameter and view voters 
as having or lacking conviction if 7 > 2 or 7 < 2, respec- 
tively. When 7 > 2, the population is quickly driven to 
consensus, while when 7 < 2, there is a bias toward the 
zero-magnetization state with equal densities of voters of 
each type. In the former case, the collective evolution 
is more extreme than the evolution of an individual. In 
the latter case, collective evolution is hindered and the 
approach to final consensus is anomalously slow. In both 
cases, we can determine the exit probability, namely, the 
probability the system ultimately reaches f consensus as 
a function of the initial composition of the population. 

In the framework of the Ising-Glauber model [?'], the 
transition rate of a voter at site i, whose states we now 
represent by (Tj = ±1, is 

w({o-}->{<J }i) = 7 + 2 - 7CT i (er,;_i + 70-j+i) 

+ (7- 2)ct < _ 1 ct <+ i, (1) 

with {a} denoting the state of all voters, {cr'}i the state 
where the i th voter has flipped, and where an irrelevant 
overall constant has been absorbed in the time scale. The 



2 



fact that the <7,<7j-|-i and the atai-i terms have the same 
coefficient is due to the left/right symmetry of the dy- 
namics. Note also that for 7 = 2, the <Ji-i<Ji + i term 
vanishes, and the equation of motion of the classic voter 
model is recovered. As shown in Glauber's original work 
and as we show below, the cji_i(Ji + i term couples the 
rate equation for the mean spin to 3-body terms. 

Because the average magnetization is not conserved, 
the model is not exactly soluble. Nevertheless, we can 
construct an approximate solution by truncating the hi- 
erarchy of rate equations of higher-order correlation func- 
tions by a simple decoupling scheme. The mean spin, 
s j = (oj) = E{ CT > o-jP{{a-};t) evolves as 



ds 



- W ({,H ff '},)F(W;4 (2) 



After straightforward steps, Eq. @ reduces to 
ds ■ 

■ = 2 7 (s j+ i + - 2(7 + 2) Sj 

-2(j-2)(a j _ l a j a j+l ), (3) 



dt 



which depends only on the mean spins Sj, Sj_ 1 and Sj+i 
when 7 = 2, as expected, but is otherwise coupled to 
higher order correlations. 

Let us first consider the mean- field limit, where the 
spins of neighboring nodes are uncorrelated. Because of 
spatially homogeneity, (sj) are all identical and we write 
the magnetization as to = (sj). Then assuming that 



{(Tj-KTjaj+i) 



m 



dm 



Eq. ^ simplifies to 



2(7- 2)(m- to 3 ), 



(4) 



which shows that the magnetization is not conserved 
when 7^2. The stable solutions of Eq. ^ are either 
consensus (to — ±1) when 7 > 2, or stasis, with equal 
densities of the two types of voters (to = 0), when 7 < 2. 

In one dimension, however, the mean-field assumption 
is not justified. To determine the behavior of the system 
when 7 7^ 2, we therefore use the approach developed 
in Q (see also to truncate the hierarchy of equa- 
tions for multi-spin correlation functions. Consider the 
rate equation for the nearest-neighbor correlation func- 



tion (CTj (7j + i) 



d(oj°>i) 
dt 



-2(7 - 2) [(cTj-io-j) + (a j+ ia j+2 ) 



+ 27 [(C7j-icr i+ i) + (cr jC r 



3 a 3+V 



+ 4 7 -4(7 + 2)(a J (7, + i). 



(5) 



To close this equation, we need to approximate the 
second-neighbor correlation function (0^0^+2). Consider 
domain walls — nearest-neighbor anti-aligned voters — 
whose density is p = (1 — ((7j<7j +1 ))/2. According to the 
transition rate in Eq. {1} , an isolated domain wall diffuses 



freely for any 7. However, when two domain walls are ad- 
jacent, they annihilate with probability P a — 7/(2 + 7) 
or they hop away from each other with probability = 
2/(2 + 7). Thus when 7 > 2, P a > and adjacent do- 
main walls have a tendency to annihilate, while they are 
repelled from each other when 7 < 2. The interaction of 
two domain walls is therefore equivalent to single-species 
annihilation, A + A — > 0, but with a reaction rate that is 
modified compared to freely diffusing reactants because 
of this interaction. Nevertheless, the domain wall density 
asymptotically decays as t~ x / 2 for any 7 < 00, and with 
an interaction-independent amplitude [9(. 

Since domain walls become widely separated at long 
times, we therefore approximate [f| the second-neighbor 

We also 



correlation function as 



(a ia 



define 1112 = ( cr i lT j+i) f° r a spatially homogeneous sys- 
tem. Now the rate equation ([5]) for the nearest-neighbor 
correlation function becomes 



<9to 2 



47 — 47TO2 



(6) 



For the uncorrelated initial condition, 777,2(0) = to(0) 2 , 
the solution is 



m 2 {t) = 1 + [to(0) 2 - 1] e 



-4yt 



(7) 



where to(0) = (sj(0)) is the average magnetization at 
t = 0. 

In a similar spirit, we decouple the 3-spin correlation 
function (crj—i 0^0^4.1) « 777777-2 [E[ and average over all 
sites, to simplify the rate equation ([3]) for the mean spin 
in a spatially homogeneous system to 



dm 
~dt 



= 2(2--f)me-^ t (m{0) 2 - 1). 



(8) 



Solving this equation and taking the t — > 00 limit, we 
obtain a non-trivial relation between the final magneti- 
zation 777(00) and the initial magnetization m(0): 



to(oo) = TO(0)e( 2 ^)(" l (°) 2 - 1 )/^ 



(9) 



It is important to realize that the average magnetiza- 
tion in a finite population does not perpetually fluctu- 
ate around this asymptotic value but rather ultimately 
reaches ±1 because consensus is the only absorbing state 
of the stochastic dynamics. We characterize this ap- 
proach to consensus by the exit probability £(x,N), de- 
fined as the probability that a population of N voters 
ultimately reaches f consensus when there are initially 
77 = xN ] voters. Since the density of f voters is 
x = (1 + to)/2 and to(oo) = 2£(x) - 1, Eq. © leads 
to 



£{x) 



(2x - l)e 2x< - 2 -^ x - 1 ^-> + I 



(10) 



The exit probability £{x) is independent of TV, but has a 
non-trivial dependence on the initial condition. Similar 
behaviors have been previously found in other opinion 
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FIG. 2: Exit probability £(x) from Eq. (flu) as a function of 
the initial density of j voters a; for different values of 7. The 
arrow indicates the position of the curves for increasing values 
of 7, namely from zero magnetization £(x) = 1/2 (7 = 0) to 
more and more consensual systems. 
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FIG. 3: Exit probability £(x) for a one dimensional system 
composed of 25 voters, for conviction parameter 7 = 0.3 and 
10. 



evolution models where the average magnetization is not 
conserved, such as the majority rule model [lol ] and the 
Sznajd model [Tlj |. 

The exit probability in (|10[) exhibits a qualitative 
change as 7 passes through 2 (see Fig. [2]). At 7 = 2, 
the pure voter model result £{x) — x, that follows from 
magnetization conservation, is recovered. When 7 < 2, 
the dynamics drives the system toward zero magnetiza- 
tion before consensus is ultimately reached. Thus the 
exit probability becomes progressively more independent 
of the initial magnetization as 7 — > 0. When 7 > 2, in 
contrast, the system is driven toward consensus; in the 
7 — ► 00 limit, £ (x) reduces to 



£(x) 



(2x-l)e" 2a;(a; - 1) + l 



(11) 



We checked the validity of Eq. (JTUJ) by simulations. To 
do so, we focused on small systems (N — 25) and di- 
rectly measured the probability £(x) that the population 



ultimately reaches a f consensus when the proportion of 
initially j voters is x by averaging over 5000 realizations 
of the dynamics. The theoretical results are in excellent 
agreement with our simulation results (see Fig. 3). 

It is in principle possible to generalize the model to 
higher dimensions d. However, in that case, each node is 
surrounded by 2d neighbors and 2d variables are there- 
fore required, i.e., the probabilities pi that an individual 
with i € [1, 2d] disagreeing neighbors changes opinion, in 
order to specify completely the voters dynamics. Such a 
general analysis, that goes in the direction of the general 
model of contagion introduced by Dodds and Watts [l2j 
and would include known models for opinion formation 
such as threshold models [l3j], will be considered in de- 
tail elsewhere. One should stress, however, that a much 
richer phenomenology can occur when the number of pa- 
rameters pi is increased [14J . In the case when each node 
is surrounded by four neighbors, for instance, and within 
a mean-field description, it is easy to show that the equa- 
tion of evolution for the average density of 1 voters is 



dx 

at 



4 

£ 

i=l 



Pi [(1 



x) 5 ~V 



f (l -*)«], (12) 



where we have again assumed that po = 0. The sta- 
ble stationary solutions of Eq. (fT2")) are either consensus 
(when pi < P4/4), a state of zero magnetization (when 
—3p\ — 3^2/2+753 + 3^4/4 < 0), or other asymmetric solu- 
tions, e.g., when pi = 1/2, pq = P2 = and ps, = p\ = 1, 
the stable solutions are x = 

Let us now return to our analysis of one-dimensional 
systems. The approach developed above can also be ap- 
plied to another simple opinion dynamics model in one 
dimension, namely, the Sznajd model [lBj]. This model is 
an appealing realization of the concept of social valida- 
tion, namely, that agents are only influenced by groups 
{e.g., pairs) of aligned voters and not by single individ- 
uals. The Sznajd is defined by the following evolution 
rule: (i) pick a pair of neighboring voters i and i + 1; (ii) 
if these voters have the same opinion it, = rjj+i, convert 
the opinion of the neighbors i — 1 and i + 2 on cither 
side of the initial pair: Oi-i = Oi = (Ji+i = Oi+2] (iii) 
repeat these steps ad infinitum or until a finite system 
necessarily reaches consensus. 

It is straightforward to show that the transition rate 
of a voter at site i is 

w{{a}^{cr }i) = -~ [cr l (cr i _2 + + + cr i+2 ) 

— crj_i(Tj_2 — o-i+io-j+2 — 2] ; (13) 

notice that the rate for changing rji involves terms that 
include £^±2- Following the same steps as in our previ- 
ous example of the voter model with variable conviction, 
the resulting equations for the average magnetization m 
and the nearest-neighbor correlation to 2 for a spatially 
homogeneous system are 



dm 
~8t 



m — mm,2 



<9m 2 



1 — rrii , 



(14) 
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FIG. 4: Exit probability £(x) for the Sznajd model. The 
system is one dimensional and composed of 25 and 100 voters 
respectively. Results are averaged over 5000 realizations of 
the random process. 



where we have again assumed that {<Jjcr 



and 



1 3° 3+1) 



-l a 3 a 3 + l) 



mrri2 ■ Here 



\°3-l a 3 U 3 + l (J 3+V 



is the correlation between the states of four contiguous 
voters. This new term is due to the distant interactions 
between voters, and is also present in the majority rule 
model. Following Q, we apply the Kirkwood approx- 
imation to factorize the 4-point function as the prod- 
uct of 2-point functions, 7714 w m 2 ,; this approach has 
proved quite useful in a variety of applications to re- 
action kinetics [16l ]. This truncation then allows us to 

solve the second of Eqs. (fL4"| to give m 2 = gnzEgf i where 
C = (m(0) 2 - l)/(m(0) 2 + 1), and finally to solve the 
equation for the average magnetization 



2m(0) 
~C l + m(0) : 



(15) 



After some straightforward algebra, we then find the exit 
probability £ (x) 



£(x) 



1 - 2x + 2a; 2 



(16) 



One should stress that such a non-trivial exit probabil- 
ity has been observed previously in simulations of the 
Sznajd model [ll|, which clearly shows that the Sznajd 
model has a tendency toward consensus and that the pre- 
diction (|16p is in remarkably good agreement with our 
simulations (see Fig. 0]). 

In this work, we investigated a general class of non- 
linear voter models of opinion dynamics in one dimen- 
sion. We considered the situation where the transition 
rate for each voter depends in a non-trivial way on the 
number of disagreeing neighbors. In general, the aver- 
age magnetization is not conserved in these models and 
the evolution of the average opinion is coupled to higher- 
order opinion correlations. It is possible to truncate the 
hierarchy of equations for these correlations in a simple 
and plausible manner. From this approach, we find that 
the system coarsens, albeit differently than in the pure 
voter model because of the interactions (repulsion or at- 
traction) of neighboring domain walls. The probability 
to reach the final state of | consensus is also shown to 
have a non-trivial initial state dependence, a feature that 
reveals the tendency of the system to inhibit or to reach 
consensus. The decoupling approximation presented in 
this letter appears to be both efficient and robust, and 
therefore it could well be useful in other models of opin- 
ion formation, language dynamics, etc. [13], where the 
interactions between agents are not conservative. 
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